Aperiodic quantum XXZ chains: Renormalization- group results 
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We report a comprehensive investigation of the low-energy properties of antiferromagnetic quantum XXZ spin 
chains with aperiodic couplings. We use an adaptation of the Ma-Dasgupta-Hu renormalization-group method 
to obtain analytical and numerical results for the low-temperature thermodynamics and the ground-state corre- 
lations of chains with couplings following several two-letter aperiodic sequences, including the quasiperiodic 
Fibonacci and other precious-mean sequences, as well as sequences inducing strong geometrical fluctuations. 
For a given aperiodic sequence, we argue that in the easy-plane anisotropy regime, intermediate between the 
XX and Heisenberg limits, the general scaling form of the thermodynamic properties is essentially given by the 
exactly-known XX behavior, providing a classification of the effects of aperiodicity on XXZ chains. We also dis- 
cuss the nature of the ground-state structures, and their comparison with the random-singlet phase, characteristic 
of random-bond chains. 

PACS numbers: 75.10.Jm, 75.50.Kj 
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I. INTRODUCTION 

At low temperatures, the interplay between lack of 
translational invariance and quantum fluctuations in low- 
dimensional strongly-correlated electron systems may induce 
novel phases with peculiar behavior. In particular, random- 
ness in quantum spin chains may lead to Griffiths phases, 1 ' 2 
random quantum paramagnetismpiii 5 . large-spin formation, 6 - 7 
and random-singlet phases. 8 ' 9 On the other hand, studies on 
the influence of deterministic but aperiodic elem ents on sim - 
ilar systems (see e.g. Refs. llMl Ill2ll3ll4ll5llr5ll7lil) . 
inspired by the experimental discovery of quasicrystals, 19 
have revealed strong effects on dynamical and thermodynamic 
properties, although much less is known concerning the pre- 
cise nature of the underlying ground-state phases^ 

Prototypical models for those studies are spin-!/2 antiferro- 
magnetic XXZ chains described by the Hamiltonian 



H = ^J,{SJSi 



(1) 



where /, > and the 5, are spin operators. In the uniform 
case (Ji = J), the ground state for chains with — 1 < A < 1 
is critical, 20 exhibiting power-law decay of the pair correla- 
tions as a function of the distance between spinsr 1 ^ as well 
as gapless elementary excitations. Such critical phase is un- 
stable towards dimerization, i.e. the introduction of alternat- 
ing couplings /odd and /even, in the presence of which a gap 
opens between the (now localized) ground state and the first 
excitated states. 22,23 - 24 This instability hints at the profound 
effects produced by fully breaking the translational symmetry 
of the system. 

Random-bond versions of these chains have been much 
studied by a real-space renormalization-group (RG) method 
introduced 25 ^ by Ma, Dasgupta and Hu (MDH) for the 
Heisenberg chain (A = 1) and more recently extended by 
Fisher, 1 - 2 - 8 who gave evidence that the method becomes 
asymptotically exact at low energies. In the last few years, the 
method has been applied and adapted to a variety of random 
systems (see e.g. Refs. I6l7l27l28l29l30l3 ll32l33l) . The ba- 
sic idea is to decimate the spin pairs coupled by the strongest 



bonds (those with the largest gaps between the local ground 
state and the first excited multiplet), forming singlets and in- 
ducing weak effective couplings between neighboring spins, 
thereby reducing the energy scale. For XXZ chains in the 
regime — 1/2 < A < 1, the method predicts the ground state 
to be a random-singlet phase, consisting of arbitrarily distant 
spins forming rare, strongly-correlated singlet pairs. 8 

Another way of breaking the translational symmetry is 
suggested by analogies with quasicrystals. These are 
structures which exhibit symmetries forbidden by tradi- 
tional crystalography, and which correspond to projections 
of higher-dimensional Bravais lattices onto low-dimensional 
subspacesi^i A one-dimensional example is provided by 
the Fibonacci quasiperiodic chain, obtained from a cut-and- 
project operation on a square lattice. 35 In this chain, spins are 
separated by two possible distances, a and b, whose sequence, 
starting from the left end of the chain, is abaab . . . This se- 
quence can be generated by repeatedly applying a substitu- 
tion (or inflation) rule a — > ab, b a, starting from a single 
distance a. Associating with each a a coupling /„ and with 
each b a coupling J/, we obtain a spin chain with couplings 
following a Fibonacci sequence. More generally, we can pos- 
tulate a two-letter substitution rule, build the corresponding 
letter sequence, and associate couplings with letters to obtain 
spin chains whose couplings follow aperiodic but determinis- 
tic sequences. 60 Quasiperiodic sequences are characterized by 
a Fourier spectrum consisting of Bragg peaks, but more com- 
plex spectra (such as singular continuous) can be generated 
by substitution rules. 36 In this work we apply the term 'ape- 
riodic' when referring to nonperiodic, self-similar sequences, 
also encompassing those which are strictly quasiperiodic in 
the above sense. 

In XX spin chains (A = 0), the low-temperature thermody- 
namic behavior can be qualitatively determined for virtually 
any aperiodic sequence by an exact RG method. 17 The ef- 
fects of aperiodicity depend on topological properties of the 
sequence. If the fraction of letters a (or b) at odd positions 
is different from that at even positions (i.e. if there is aver- 
age dimerization), then a finite gap opens between the global 
ground state and the first excited states, and the chain becomes 
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noncritical. Otherwise, the scaling of the lowest gaps can be 
classified according to the wandering exponent co measuring 
the geometric fluctuations g related to nonoverlapping pairs 
of letters, 17 which vary with the system size N as g ~ N w . If 
CO < 0, aperiodicity has no effect on the long-distance, low- 
temperature properties, and the system behaves as in the uni- 
form case, with a finite susceptibility at T — 0. If co = 0, as 
in the Fibonacci sequence, aperiodicity is marginal and may 
lead to nonuniversal power-law scaling behavior of thermo- 
dynamic properties. If CO > 0, aperiodicity is relevant in the 
RG sense, affecting the T = critical behavior and leading 
to exponential scaling of the lowest gaps A at long distances 
r, according to the form A ~ exp(— r a ). In particular, for 
sequences with CO = 1/2, geometric fluctuations mimic those 
induced by randomness, and the scaling behavior is similar to 
the one characterizing the random-singlet phase. 8 

In contrast, results for the effects of aperiodicity on low- 
energy properties of XXZ chains have been so far scarce, 
and restricted to particular sequences. Vidal, Mouhanna and 
Giamarchi 13,14 studied the related problem of an interact- 
ing spinless fermion chain with Fibonacci or precious-mean 
potential by using bosonization techniques, which are valid 
in the weak-modulation regime (J a ~ J/,). At half filling, 
where the system corresponds to an XXZ chain in zero ex- 
ternal field, their calculations predict that aperiodicity will 
drive the system away from the usual Luttinger-liquid behav- 
ior for < A < 1 . A similar conclusion is drawn from stud- 
ies on a Hubbard chain with hoppings following a Fibonacci 
sequence. 16 Density-matrix renormalization-group (DMRG) 
results on XXZ chains with precious-mean couplings^ and 
recent real-space RG calculations on the Fibonacci XXZ 
chair. 37 - 38 (also based on the MDH scheme), likewise pre- 
dict that low-temperature properties are different than in the 
uniform chains. The zero-temperature magnetization curve of 
Fibonacci XXZ chains has also been investigated, 18 with em- 
phasis on determining the plateau structure. 

Our aim in this paper is to investigate the effects of arbi- 
trary aperiodic coupling distributions on the low-temperature 
properties of XXZ chains, reinforcing and extending our previ- 
ous results. 39 From an adaptation of the Ma-Dasgupta-Hu RG 
scheme, we obtain information about low-temperature ther- 
modynamics and ground-state correlation functions for sev- 
eral aperiodic sequences. Our results, which are presumably 
exact in the strong-modulation limit (J a <C Jb or J a ^> Jb), 
point to the following conclusions: 

• The exact classification found in the XX limit can ar- 
guably be extended to XXZ chains in the anisotropy 
regime < A < 1 . We predict that dimerized aperiodic- 
ity opens a gap to the lowest excitations, and that other- 
wise the effects of aperiodicity on the low-temperature 
thermodynamics are gauged by the same exponent co, ir- 
respective of anisotropy. In particular, sequences which 
are strictly marginal in the XX limit continue to be so for 
anisotropics < A < 1, but may be marginally relevant 
in the Heisenberg limit. 

• On the other hand, co is found not to define the behavior 
of correlation functions, although ground-state struc- 



tures in the presence of marginal or relevant couplings 
also reflect self-similar properties of the sequences. 
Dominant correlations correspond to well-defined dis- 
tances, related to the rescaling factor of the sequences 
(contrary to the random-singlet phase, where no such 
characteristic distances exist), and two types of be- 
havior are possible: either the chains can be decom- 
posed into a hierarchy of singlets, forming a kind of 
'aperiodic-singlet phase', or into a hierarchy of effec- 
tive spins, in which case low-energy excitations involve 
an exponentially large number of spins. This is in sharp 
contrast both to the gapless spin-wave excitations of the 
uniform chains and to the gapped triplet-wave excita- 
tions of the dimerized chains. 

• Based on second-order calculations, the long-distance 
decay exponents of average ground-state correlation 
functions are seen to vary with the coupling ratio in 
the presence of strictly marginal aperiodicity. Other- 
wise, strong universality (i.e. independence of the expo- 
nents on both coupling ratio and anisotropy) is obtained 
for the whole line < A < 1, although different decay 
exponents may emerge in the Heisenberg limit. Also, 
the scaling form of typical (rather than average) cor- 
relations follows essentially the same scaling form as 
the energy gaps, similarly to what happens for random- 
bond chains. 

In order to make the paper self-contained, we begin by re- 
viewing some known results. So, in Sec. |II]we present the 
basics of the Ma-Dasgupta-Hu scheme, as applied to random- 
bond XXZ chains, and summarize the properties of the un- 
derlying random-singlet phase. Also, in Sec. |III]we provide a 
short discussion on aperiodic sequences, as well as a sketch of 
the exact RG results for XX chains with aperiodic couplings. 
Our adaptation of the Ma-Dasgupta-Hu method to aperiodic 
XXZ chains is described in Sec. II VI and results for marginal 
and relevant aperiodicity are presented in Sees. IVland lVTl The 
final section is devoted to a discussion and conclusions. There 
are also two appendices, in which some important technical 
points are detailed. 



II. RANDOM-BOND SPIN CHAINS AND THE 
MA-DASGUPTA-HU METHOD 

Consider an antiferromagnetic quantum spin-!/2 chain de- 
scribed by the Hamiltonian 

H = £j,(S*S* +l +S>Sl l+ A,Stf +l ), (2) 

i 

where 7, > and all anisotropics are such that < A,- < 1. 
Let us assume that the couplings /, are randomly distributed 
according to a broad probability distribution having an 
upper cutoff 7 max . Under such conditions, in a finite but large 
chain, there is a strongest bond Jo ~ 7 max connecting, say, 
spins S\ and S2, which on their turn are coupled to their other 
nearest neighbors Si and S r by weaker bonds J\ and J r . The 
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local Hamiltonian connecting S\ and S2 is 



III. APERIODIC SEQUENCES AND XX CHAINS 



#0 = Jo (S1S2 



- S\S' 2 + AqSjSj) 



whose ground state is a singlet, separated from the first ex- 
cited states by an energy gap Ao = \{l +Ao)7o- The idea 
behind the Ma-Dasgupta-Hu2^2^ method is that, at tempera- 
tures below Ao, Si and S2 can be decimated out of the system, 
since they couple into a singlet, giving a negligible contribu- 
tion to thermodynamic properties. Nevertheless, their virtual 
excitations induce a weak effective coupling between Si and 
S r , described by the Hamiltonian 



H' 



■f {s;s; 



The parameters J' and A' can be obtained by second-order per- 
turbation theory (see AppendixlAl. and are given by 



1 



1+A Jq 



and A' 



1+A 



A/A r 



(3) 



Notice that /' is smaller than either //, J r or Jo; likewise, un- 
less all Af = 1, A' is smaller than either A/ or A r . Thus, af- 
ter eliminating Si and S2, the overall energy scale is reduced. 
The previous steps can be repeated with the next largest bond, 
which most probably is not J'. 

Starting from Eqs. Q, Fisher 8 was able to write and solve 
recursion relations for the probability distribution of the ef- 
fective couplings. The fixed-point distribution is presumably 
independent of the initial couplings, 40 and diverges as a power 
law for /' — > + , indicating that the perturbative approach 
leading to Eqs. Q becomes essentially exact for asymptot- 
ically low energies. 

The ground state is a 'random-singlet' phase, consisting of 
arbitrarily distant spins forming rare, strongly correlated sin- 
glet pairsi^i Exciting a singlet whose spins are separated by 
a distance r costs an energy of order A, with a dynamic scaling 
form 



where p and ro are constants. At low temperatures, the zero- 
field susceptibility diverges as 



1 



rin 2 r' 



Average ground-state correlations are dominated by the rare 
singlet pairs, and decay as a power law, 



C(r) 



1 



where the bar denotes average over the whole chain, while 
typical correlations are short-ranged, following 



CtypM 



These asymptotic results are independent of the anisotropies 
A,, as long as < A; < 1 for all ;, with the same distribution 
on even and odd bonds. 



Following closely the analysis of Hermisson, 17 in this Sec- 
tion we consider antiferromagnetic quantum XX chains, de- 
scribed by the Hamiltonian 



(4) 



where now the strengths of the site-dependent couplings 7, can 
be either J a or Jj,, and are distributed according to determin- 
istic but aperiodic binary sequences, obtained by substitution 
(or inflation) rules of the form 



o : 



w a 

Wb 



where w a and Wb are words (finite strings) composed of letters 
a and b. A well-know example is provided by the Fibonacci 
sequence, whose substitution rule is 



■ ab 
a 



Starting from a single letter a, repeated application of Oft, 
yields strings with lengths given by the Fibonacci numbers 
1, 2, 3, 5, 8,..., ultimately producing a letter sequence 
abaababaaba . . ., for which no period can be identified. 

Given an inflation rule o, various statistical properties 4 ^ of 
the associated sequence are enclosed in the substitution matrix 



M: 



#a{w a ) #a{wb) 
#b(Wa) #b{wb) 



where # a {wb) denotes the number of letters a in the word Wb- 
The largest eigenvalue of M, X+, gives the asymptotic scaling 
factor of the string length (i.e. the ratio between the lengths of 
the strings corresponding to successive iterations of the rule 
a); the entries of the corresponding eigenvector are propor- 
tional to the frequencies p a and pb of letters a and b in the 
limit (infinite) sequence. 

The remaining eigenvalue, A,_, is related to the geometric 
fluctuations of the sequence, which are defined in the follow- 
ing way. Let A^' be the number of letters a in the string ob- 
tained after n iterations of o, and N n be the corresponding total 
number of letters (the length of the string). Then, a measure of 
the geometric fluctuations induced by the sequence is the dif- 
ference g n between N" and the number of letters a expected 
from the limit-sequence frequency p a , and this behaves as 

gn = \K- Pa N n \~\X-\ n . 

Since ~ A/J_, this last equation can be rewritten as 

8n *~ A^ n ' , 
by defining the 'wandering exponent' 

ln|A,_| 



CO/ 



In Am 
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If CO/ < 0, fluctuations become smaller as the string grows, 
and the sequence looks more and more 'periodic'. On the 
other hand, if CO/ > 0, fluctuations increase without limit. 
The marginal case CO/ = is in general connected to loga- 
rithmic fluctuations. It can be shown 43 that substitutions for 
which CO/ < generate quasiperiodic (or limit-quasiperiodic) 
sequences. 

The concept of geometric fluctuations is essential in estab- 
lishing the Harris-Luck criterion^ii 4 ^ for the relevance of in- 
homogeneities to the critical behavior of magnetic systems. 
The criterion states that, if fluctuations in the local parameters 
controlling the criticality of the system vary with some char- 
acteristic length L as g ~ L m , then there is a critical value of 
the exponent CO above which the presence of inhomogeneities 
can affect the critical behavior. This happens for 



co > 1 



1 



(5) 



where d is the number of dimensions along which inhomo- 
geneities are distributed, and v is the correlation-length crit- 
ical exponent of the underlying uniform system. Randomly 
distributed inhomogeneities lead to CO = 1/2, and the general 
Harris criterion 44 ' 46 ' 47 is recovered. 

The ground-state of the model in Eq. @ is critical in the 
uniform limit (7,- = J): there is no energy gap to the low-lying 
excitations, and pair correlations decay as power laws^ 



C°"(r) = (SfS^ 



with r) xx = rp = 1/2 and r|" = 2. This phase is unstable to- 
wards dimerization (i.e. the presence of couplings J and J e 
alternating between odd and even bonds), in which case a gap 
opens in the low-energy spectrum, and ground-state correla- 
tions become short-ranged. More generally, the model ex- 
hibits a (zero-temperature) quantum phase transition between 
two dimer phases forii 4 ^ 



8 = ln72j-i — ln/2/ = 0. 



(6) 



If /2/-I = and J%j = J e , the phase transition occurs for 8 = 
ln(7 /7 e ) = 0, and belongs to the Onsager universality class, 
with v = 1 . 

When the couplings Jj are chosen according to aperiodic 
sequences for which the fractions of letters a (or b) at even 
and odd positions are different, Eq. is not satisfied, and 
the system is in a dimer phase. This suggests that the local 
parameters defining the criticality of XX chains are the shifts 
8 ; =ln(72/-i/^2/)- In order to study the fluctuations of the 8 ; -, 
which depend on two consecutive couplings, we must usually 
consider the sequence of nonoverlapping letter pairs associ- 
ated with a given aperiodic sequence. To build the inflation 
rule c^ 2 ' for such pairs, it is necessary to iterate the original 
rule a until the strings obtained from a single a and b have 
lengths of the same parity. As an illustration, let us take the 
Fibonacci sequence. Applying Oft, three times yields 



J fb • 



abaab 
aba 



Noting that the pair bb does not occur in the sequence, we 
readily obtain 



(ab) (aa){ba) (ba) (ab) 

(ab)(aa)(ba)(ba) 

(ab)(aa)(ba)(ab) 



For a general pair inflation rule a' 2 ', we can define an asso- 
ciated substitution matrix 




M 



(2) _ 



/ #aa{Waa) #aa(w a b) #aa(w b a) #aa(wbb) 

#ab(w aa ) #ab(w a b) #ab(w b a) #ab(wbb) 

#ba(w aa ) #ba(w a b) #ba(wba) #ba(wbb) 

V #bb{w aa ) #bb{w a b) #bb(w b a) #bb(w b b) 



where now # a b(Wba) denotes the number of pairs ab in the 
word associated with the pair ba. The leading eigenvalues A-i 
and A-2 of M^ 2 ' define another wandering exponent 



co = 



\n\X 2 
lnA-i 



which governs the fluctuations of the letter pairs, and conse- 
quently of the 8j. It is essential to note that CO is in general 
different from CO/: for the Fibonacci sequence, for instance, 
we have CO/ = —1, but co = 0& 

By an exact renormalization-group treatment, Hermisson 17 
was able to build recursion relations for effective couplings, 
and to show that, in agreement with the above heuristic argu- 
ment, the eigenvalues A., of M' 2 ' give directly the RG eigen- 
values yi around the uniform fixed point of XY chains with 
aperiodic couplings, 



yi 



\n\Xj\ 
InA-i 



while the corresponding eigenvectors yield the scaling fields. 
Thus, aperiodicity is relevant in the RG sense (i.e. it moves 
the RG flows away from the uniform fixed point) if the next- 
to-leading eigenvalue y 2 = CO is positive, exactly as predicted 
by the Harris-Luck criterion, Eq. 10, with d = v = 1 

For a large class of aperiodic sequences fulfilling Eq. 0, 
X2 is given in the XX limit by an integer k. When a pair sub- 
stitution rule can be defined, this integer is simply given by 

k = #ab{Wab) - #ba(w a b)- 

Thus, the wandering exponent in the XX limit is of the form 

ln/t 



co = 



lnT' 



(7) 



where x = X\ corresponds to the rescaling factor of the se- 
quence of letter pairs. 

It is also possible to determine the scaling of the lowest en- 
ergy levels A with the system size r. For irrelevant or marginal 
aperiodicity (co < 0) we have 

A~r-\ (8) 

with a dynamical exponent z equal to unity if CO < 0, but which 
can vary with the coupling ratio J a /Jb m the marginal cases 
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(co = 0). Relevant aperiodicity (co > 0) leads to a different 
scaling form, 

A ~ exp (— /jr m ) , (jj = constant) (9) 

and to a formally infinite dynamical exponent. From Eqs. (|8j 
and scaling forms for low-temperature thermodynamic 
properties such as the specific heat and zero-field susceptibil- 
ity can be obtained, as discussed in the next sections. Ground- 
state correlation functions, however, do not seem to be simply 
accessible from the exact RG treatment. 

Since the critical phase of uniform XXZ chains is also un- 
stable towards dimerization in the whole anisotropy regime 
< A < 1, one might expect that the relevant geometrical 
fluctuations in the presence of aperiodic couplings would be 
somehow related to the 8 ; - defined above. 63 Consequently, the 
exponent co would be involved in determining the scaling be- 
havior of thermodynamic properties of aperiodic chains for 
all anisotropies intermediate between the XX and Heisenberg 
limits. The results of the next sections indeed provide evi- 
dence that this seems to be the case. 



IV. THE MA-DASGUPTA-HU METHOD FOR APERIODIC 
XXZ CHAINS 

We now wish to investigate the effects of aperiodic cou- 
plings on XXZ chains described by the Hamiltonian in Eq. 0. 
Based on the success of the Ma-Dasgupta-Hu scheme in pre- 
dicting the properties of random-bond chains, we expect that 
it also works in the presence of aperiodicity. We concentrate 
on the case of uniform anisotropy (A; = A), but more general 
situations can be considered. 

Applying the MDH method to aperiodic chains requires 
taking into account that now, since we have only two dis- 
tinct coupling constants, there are many spin blocks with the 
same (largest) gap at a given energy scale. Also, those blocks 
may consist of more than two spins, in which case effective 
spins would form upon renormalization. The strategy is to 
sweep through the lattice until all blocks with the same gap 
have been renormalized, leading to new effective couplings 
(and possibly spins). Then we search for the next largest gap, 
which again corresponds to many blocks. When all possible 
original blocks have been considered, there remains some un- 
renormalized spins, possibly along with effective ones, defin- 
ing new blocks which form a second generation of the lattice. 
The process is then iterated, leading to the renormalization of 
the spatial distribution of effective blocks (or bonds) along the 
generations. 

Due to the self-similarity inherent to aperiodic sequences 
generated by inflation rules, it is natural that the block dis- 
tribution reaches a periodic attractor (usually a fixed point or 
a two-cycle) after a few lattice sweeps; numerical implemen- 
tations of the method indicate that this attractor is indepen- 
dent of the anisotropy A for all coupling ratios. By study- 
ing recursion relations for the effective couplings, we can ob- 
tain analytical results. As the RG steps proceed, the cou- 
pling ratio usually gets smaller, suggesting that the method 
becomes asymptotically exact. This picture holds for marginal 



(co = 0) and relevant (co > 0) aperiodicity. Irrelevant ape- 
riodicity is characterized by a wandering exponent CO < 0, 
meaning that geometric fluctuations become negligible at long 
distances. An example is provided by the Thue-Morse se- 
quence, generated by the substitution rule a — ► ab, b — > ba, 
for which CO = — °°. Applying the MDH scheme to the Thue- 
Morse sequence leads to an effective coupling ratio which ap- 
proaches unity along the generations, although the couplings 
themselves become smaller. This means that the perturbative 
approach in the core of the MDH scheme eventually breaks 
down, and no asymptotic behavior can be obtained. However, 
this intuitively agrees with the picture that irrelevant aperiod- 
icity leads to the same critical properties as the uniform model, 
where all couplings have the same value. 

For sequences where the fraction of letters a (or b) at odd 
bonds is different from that at even bonds (i.e. where the se- 
quence induces average dimerization), one generally expects 
that a finite gap opens between the global ground state and the 
first excited states, independent of the value of the wandering 
exponent CO. This is the case of the period-doubling sequence, 
built from the substitution rule a — > ab, b — > aa. Upon appli- 
cation of the MDH method, after a few lattice sweeps (with 
the precise number depending on the strength of dimeriza- 
tion) we reach a situation where, say, all strong bonds occupy 
even positions, whereas all bonds at odd positions are weaker. 
Thus, all remaining couplings are necessarily decimated in a 
last lattice sweep, generating a final effective coupling which 
approaches zero exponentially with the system size N. This 
can be interpreted as indicating that there is no correlation be- 
tween spins separated by large distances, in agreement to what 
happens in gapped Heisenberg and XX chains. In the pres- 
ence of average dimerization, few quantitative predictions can 
be drawn from the MDH method; one of them is an estimate 
of the excitation gap, whose order of magnitude is provided 
by the value of the strong bonds in the final lattice sweep. In 
contrast, randomly dimerized Heisenberg chains in the strong- 
randomness limit are in a gapless Griffiths phase, exhibiting 
short-range correlations but a diverging susceptibility i2i 

In a general situation, the blocks to be renormalized consist 
of n spins connected by equal bonds Jo, and coupled to the rest 
of the chain through weaker bonds // and J r . As discussed in 
Appendix lAl the ground state for blocks with an even num- 
ber of spins is a singlet (as in the original MDH method), and 
at low energies we can eliminate the whole block, along with 
J i and J r , leaving an effective antiferromagnetic bond J' cou- 
pling the two spins closer to the block and given by second- 
order perturbation theory as 

, J\J r 
J =Y«-— (neven), 
Jo 

with A-dependent coefficients y„ . On the other hand, a block 
with an odd number of spins has a doublet as its ground state; 
at low energies, it can be replaced by an effective spin con- 
nected to its nearest neighbors by antiferromagnetic effective 
bonds 

J' l r = Y„7/, r (« odd), 
whose values are calculated by first-order perturbation theory. 
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In general, the anisotropy parameters are also renormal- 
ized and become site-dependent; for n even, the effec- 
tive anisotropy is A' = 8„(Ao)A/A r , while for n odd A' lr = 
8„(Ao)A /)r , with |8„(A)| < 1 forO< A< 1 andS„(l) = l. So, 
for < A < 1 the A,- flow to the XX fixed point (all A, = 0), ulti- 
mately reproducing the corresponding scaling behavior, while 
for the Heisenberg chain all A,- remain equal to unity. An ana- 
lytical treatment of the intermediate anisotropy regime is pos- 
sible (see Ref. l38h . leading to a prediction of the effective 
coupling ratio for which the system crosses over to the XX 
behavior. However, for simplicity, we present analytical cal- 
culations for the XX and Heisenberg limits, showing some nu- 
merical results for the general case < A < 1 . If we start with 
a uniform anisotropy A > 1, the A, grow without limit, and 
the system ultimately behaves like an antiferromagnetic Ising 
chain, suppressing all quantum fluctuations. For A < 0, the 
J„ coefficients for n even become larger than unity, so that, if 
the modulation is not strong enough, the MDH scheme may 
produce effective couplings which are larger than the original 
couplings, leading to 'bad' decimations; moreover, the two- 
spin local gap closes as A — > — 1 . This puts the MDH results 
under suspicion, requiring a more careful analysis which is 
beyond the scope of the present work. 

Correlation functions can be calculated at zeroth order by 
assuming that only spins which eventually appear in the same 
renormalized block are correlated. Note that an effective spin 
represents all spins in an original block via Clebsch-Gordan 
coefficients (see Appendix lAl. and this allows us to calculate 
correlations between any two spins whose effective spins end 
up in the same block at some stage of the RG process. In or- 
der to estimate correlations between other spin pairs we must 
expand the local ground states up to second order in 7/ r /7o- 
This requires lengthy calculations (see AppendixlEl. and we 
restrict applications of this expansion to the simplest yet il- 
lustrative cases of sequences where only two-spin blocks are 
involved in the RG steps. 

In the next two sections, we present a detailed discussion 
of the results obtained by applying the MDH scheme to se- 
quences inducing marginal or relevant aperiodicity. 

We should mention that similar strong-modulation pertur- 
bative approaches have been applied to investigate the spectral 
properties of noninteracting electrons with aperiodic hopping 
parameters or single-site potentials (see e.g. Refs. 15(15 ll 
However, the XXZ chain with nonzero anisotropy A is mapped 
by the Jordan-Wigner transformation onto a half-filled inter- 
acting electron system, for which, to the best of our knowl- 
edge, no such studies exist. 

V. MARGINAL APERIODICITY 

A. The Fibonacci sequence 

First we apply the method to chains with Fibonacci cou- 
plings. This is the simplest example of the quasiperiodic 
precious-mean sequences with marginal fluctuations. 17 A few 
bonds closer to the left end of the original chain, along with 
induced effective couplings, are shown in Fig. ^for J a < Jt,. 



Figure 1 : Left end of the Fibonacci XXZ chain with J a < J^. Dashed 
(solid) lines represent weak (strong) bonds, while circles indicate the 
positions of the spins. Apart from a few bonds close to the chain 
ends, the effective couplings also form a Fibonacci sequence. 

In this case, only singlets are formed by the RG process; apart 
from a few bonds close to the chain ends, the renormalized 
lattice is again a Fibonacci chain. An effective coupling J' b 
is induced between spins separated by only one singlet pair, 
while J' a connects spins separated by two singlet pairs, and in 
terms of the original couplings we have 




The bare coupling ratio is p = J a /Jb, its renormalized value 
being p' = J2P- In each generation j, all decimated blocks 
have the same size ry and gap A/ (proportional to the effective 
Jb bonds). The recursion relations for p and A are given by 

Pj+i=Y2Pj and Ay+i = Y 2 pyA;. (10) 

For the XX chain 72 = 1, and thus p ; +i = pj, corresponding 
to a line of fixed points. On the other hand, for the Heisen- 
berg chain 72 = l /2, so that pj + \ < pj, leading to a stable fixed 
point poo = 0; since the perturbative approach on which the 
MDH scheme is based works for p< 1, the method can be 
expected to yield asymptotically exact results for the Heisen- 
berg Fibonacci chains. In both cases, solving Eqs. dlOt gives 
the gap in the /'th generation in terms of the original coupling 
ratio p and gap A, 

Ay = Yfp 2; 'A. 

(Notice that, since 72 depends on the anisotropy, this last equa- 
tion is valid only for A,- = or A, = 1; in the intermediate 
anisotropy regime, the variation of 72 along the generations 
must be taken into account. 38 ) The distance between spins 
forming a singlet in the y'th generation defines a characteristic 
length rj, corresponding to the Fibonacci numbers r j = 1, 3, 
13, 55, . . .; for ;> 1 the ratio rj + i/rj approaches (|) 3 , where 
(j) = (1 + V5) /2 is the golden mean. So we have r ; - ~ ro<() 3j , 
where ro is a constant, and we obtain the dynamical scaling 
relation 

A, ~ r^e-^^ilro) 

J J 

with £ = — 2/3lnp/ln(|) and /j = -ln72/91n 2 (j). For the 
Heisenberg chain (72 = l /2), Eq. il lb describes a weakly ex- 
ponential scaling (with a formally infinite dynamical expo- 
nent), but not of the form A ~ exp(— r m ) found for the XX 
chain with relevant aperiodicity (co > 0). For the XX chain 
(72 = 1)) M = and we can identify ^ with a dynamical ex- 
ponent z, whose value depends on the coupling ratio, lead- 
ing to nonuniversal scaling behavior, characteristic of strictly 



7 



Figure 2: Left end of the Fibonacci chains with J a > J/,. Effective 
spins form in the first lattice sweep, giving rise to a Fibonacci chain 
with the roles of the weak and strong bonds interchanged, exactly as 
in the original lattice in Fig. Q 

marginal operators. (We can check that z — C corresponds to 
the asymptotic form of the exact XX expression for p <C 
1.) This nonuniversality should hold in the anisotropy regime 
< A < 1 with a 'bare' value of p defined at a crossover scale. 
Note that, taking into account the scaling form A ~ exp (— r m ) 
valid for relevant aperiodicity, we can view the above Heisen- 
berg scaling form (fj ^ 0) as a marginally relevant (co — > + ) 
case. The result in Eq. M \\ has also been obtained in Ref. |37l 
If we choose J a > Jy, blocks with three spins connected 
by two strong bonds appear in the chain, producing effective 
spins upon renormalization. However, as illustrated in Fig. 
|2j the first lattice sweep yields again a Fibonacci chain with 
the roles of weak and strong bonds interchanged, exactly as in 
Fig. ^ The effective couplings in the second generation are 
given by 

jZ 

J'a=l2%T and f b =$J b , (12) 
and the coupling ratio is now 

p= l = ^r (T2p) ' 

which is larger than one, showing that J' a < J',, Thus, we can 
apply the same analysis as in the case with J a < J/,, but now 
with 'bare' couplings given by Eq. d!2i . So, in the XX chain, 
since fi = 1, the MDH method predicts scaling forms which 
are symmetric under p — > 1/p, in agreement with the exact 
treatment. Uiii 

The susceptibility %{T) can be estimated 8 by assuming that, 
at energy scale Aj ~ T, singlet pairs are effectively frozen, 
while unrenormalized spins are essentially free, contributing 
Curie terms to the susceptibility. Thus, if nj ~ rj 1 is the num- 
ber of surviving spins in the jth generation, %(T ~ Aj) ~ 
tij + \/Aj. This already gives reasonable results, as indicated 
by comparison with those obtained for the XX chain from nu- 
merical diagonalization of finite chains, 39 based on the free- 
fermion method. 48 However, a more useful approximation can 
be obtained by noting that, in the jth generation, we can view 
the resulting lattice as composed of 'independent' singlets in 
which a pair of spins is coupled via an XXZ interaction with 
effective bond and anisotropy parameters J^jp and A, . Since 
the fraction of such singlets with respect to the number of orig- 
inal bonds is (ni — n 7+ i)/2, the free energy per site of the 
whole system, in the presence of an external field h — > 0, can 
be estimated as 

= l£^±iF pail (jj?>,4P;h,T) , (13) 




Figure 3: (Color online) Specific heat of the Fibonacci XXZ chains 
for Ja/Jb = V 10 an d mree different values of the uniform anisotropy 
A, as given by the 'independent-singlet' approximation. 



where F pa i r (/, A;h, T) is the free energy of a pair of spins in- 
teracting via the Hamiltonian 

Hpair = J (SiSf + S\S y 2 + AS\S Z 2 ) - h{S\ + S z 2 ). 

Iterating the recursion relations for the effective couplings J a j 
and A a b, we can determine their values in each generation, 
and evaluate numerically the sum in Eq. d!3> to obtain the 
free energy. Thermodynamic properties such as the zero-field 
susceptibility % and the specific heat c can be obtained by the 
relations 

a 2 / 3 2 / 

As an example, Fig. [3]shows plots of the specific heat of Fi- 
bonacci XXZ chains with J a /Jb = V 10 an d th ree values of the 
anisotropy A = A a — Aj,, corresponding to the XX and Heisen- 
berg limits and to an intermediate case (A = 9 /io). The results 
for the XX limit agree very well with those obtained from 
numerical diagonalization, although the agreement becomes 
worse for larger coupling ratios; in particular, the specific-heat 
scaling lawiiiii 

^>~^(^)- 

with z = C(p) an d G c a function of unit period, is fully sat- 
isfied, reflecting the strictly marginal character of the aperi- 
odic perturbations. This is not the case in the Heisenberg 
limit, and the logarithmic amplitudes of the oscillations in 
the specific heat become larger with decreasing temperatures, 
reflecting the weakly exponential dynamical scaling in Eq. 
(II Q . For intermediate anisotropics, there is a crossover from 
Heisenberg-like to XX-like behavior as the temperature is low- 
ered; the larger amplitude of the low-temperature oscillations 
corresponds to those of an XX chain with a 'bare' coupling 
ratio p e fl- < p defined at a crossover scale in which effective 
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Figure 5: Left end of the silver-mean chain with J a < J/,. Effective 
couplings j' b correspond to the original J a bonds, while j' a connects 
spins separated by one singlet pair. Apart from the leftmost bond, the 
effective couplings also form a silver-mean sequence. 
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Figure 4: (Color online) Ground-state correlations as a function of 
the distance between spins for the Fibonacci XX chain. The curves 
are obtained from numerical diagonalization of closed chains with 
2584 sites. For J a /Jb — V 1() Q ower curves), dominant correlations 
correspond to distances rj = 1, 3, 13, 55 and 233, for which C xx and 
C zz are nearly equal, as predicted by the MDH method (circles), and 
decay as 1 jr (dotted curve). Larger coupling ratios lead to a slower 
decay of C xx (and a faster decay of C zz ), as seen for J a /Jb = V 2 
(upper curve, offset for clarity). 



anisotropics become negligible. (For a detailed analysis, see 
Ref.lM) 

As all singlets formed in the jth generation have length rj 
and the bond distribution is fixed, the average ground-state 
correlation between spins separated by a distance rj can be 
estimated as 



Figure 6: Left end of the silver-mean chain with J a > J/,. Effective 
spins form in the first lattice sweep, producing the same structure as 
the third generation for J a < Jb- 



B. The silver-mean sequence 

The silver-mean sequence is obtained from the substitution 
rule a — > aab, b — > a, and the rescaling factor predicted in the 
XX limit isiZ 

X S m = 1 + V2- 

When the MDH scheme is applied, the first lattice sweep also 
generates a silver-mean sequence, identical to the original one 
for J a < Jh, but with the roles of weak and strong bonds in- 
terchanged for J a > Jb, as shown in Figs. [3] and [6] In the 
latter case, the second-generation structure is identical to the 
third-generation lattice obtained for J a < Jb, a situation we 
can assume without loss of generality. So, we can write the 
recursion relations 



y 2 ^- and f b =J a , 
Jb 



C™(rj)^( S ?Si 



co\(rij— =o\co\r, 1 , (14) from which we get 



where the bar denotes average over all possible pairs, o is 
a constant, a = x,y,z, and co is the correlation between the 
two spins in a singlet, given by co = — '/ 4 f° r tne Heisenberg 
chain and for both a = x and a = z in the XX chain. We point 
out that these should be the dominant correlations, and spins 
separated by distances other than r ; are predicted to be only 
weakly correlated. As shown in Fig. |4] results from numeri- 
cal diagonalization for the XX Fibonacci chain with p = l/io 
agree very well with the MDH predictions. Note that correla- 
tions in the uniform XX chair. 48 decay as C xx (r) ~ r~ 1//2 and 
C zz {r) ~ r~ 2 , so that dominant xx (zz) correlations in the Fi- 
bonacci chain are weaker (stronger) than in the uniform chain. 
Due to the strictly marginal character of the fluctuations in- 
duced by the aperiodic couplings, deviations from the predic- 
tions in Eq. H41 appear in the XX chain for larger values of 
p, as also shown in the figure. This point will be further dis- 
cussed in the next subsection, but these deviations should not 
be present in the Fibonacci Heisenberg chain, where aperiod- 
icity can be viewed as marginally relevant. 



fa 



Ja 



A' JL 



J2f=j2p and -r = -r=P- 
Jb A J h 



These are similar to the relations found for the Fibonacci 
chains. The length of singlets formed in the jth generation 
is rj = 1, 1, 3, 7, 17, 41, whose asymptotic ratio is 
r i+i / r i = x sm- Thus, solving the recursion relations yields 



-Up) e -fihx 2 (rj/r ) 



with 



C(P) 



In 



In Xsm 



and 



m(yf) 

In 2 Xsm 



(15) 



so that in the XX limit the scaling again corresponds to a 
nonuniversal power-law behavior with a dynamical exponent 

z = «p). 

As in the Fibonacci chains, pair correlations in the ground 
state can be estimated by noting that only singlets are pro- 
duced by the RG process, and we conclude that for p <C 1 the 
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Figure 7: (Color online) Ground-state pair correlations in the silver- 
mean XX (a) and Heisenberg (b) chains, for p = 1/4 (upper curves) 
and p = 1/ 10 (lower curves), obtained from the second-order MDH 
scheme (chains with 81 19 sites). Solid and dashed curves correspond 
to average and typical correlations, respectively. 



dominant correlations (those between spins separated by the 
characteristic distances rj = 1, 3, 7, 17,. . .) should behave as 



C(ry). 



1 



while correlations between spins separated by other distances 
should be negligible. 

However, for XX chains, this is a rough approximation if the 
coupling ratio is not too small, and free-fermion calculations 
reveal a power-law decay of both C xx (rj) and C zz (r,) with p- 
dependent exponents, as in Fig. |4] This can be accounted for 
by the MDH method if we expand the ground-state vector to 
second-order in p, as described in Appendix IH1 Results of 
such calculations are shown in Fig. Qfor p = [ /4 and p = '/io, 
in the XX and Heisenberg limits. Both average and typical 
correlations are plotted; the latter, defined by 



C? y «(r)=exp(ln 



(SfSf +r )\) 



filter out the contribution of those pairs of spins most strongly 
correlated, yielding an estimate of the correlation between 
two arbitrary spins separated by a distance r. In the random- 
singlet phase, characteristic of random-bond chains, 8 aver- 
age correlations decay algebraically as C(r) ~ 1 /r 2 , whereas 
typical correlations are short-ranged, following C typ (r) ~ 
exp(— y/r/ro). This is due to the fact that average correla- 
tions are dominated by the rare singlet pairs, while the corre- 
lation between a typical pair of spins is of the order of some 
intermediate effective coupling (see AppendixlEV 

As shown in Fig. EI a), this picture does not hold for silver- 
mean XX spin chains. As the coupling ratio is lowered, av- 
erage and typical correlations exhibit clearly distinct behav- 
ior, but both C xx (rj) and Cfy p (r) still follow approximately a 
power law, with p-dependent exponents, reproducing the re- 
sults of the free-fermion calculations. This nonuniversality is 
related to the marginal character of the precious-mean fluc- 
tuations, which keeps the effective coupling ratio unchanged 



along the RG process, and can be qualitatively understood 
from the following argument. For each singlet pair coupled 
by a strong bond and whose spins are separated by a charac- 
teristic distance ry, there exists a certain number of other spin 
pairs separated by the same distance ry, but connected through 
weaker bonds, whose correlation (see Appendix[5J is smaller 
than the strongest ones by factors of order p, p 2 , p 3 , etc. The 
average correlation can be estimated as 



C xx (n) ~ (l+aiP + fl2p 2 H \-ajP j )rj 



(16) 



where the a„'s are proportional to the fractions of pairs giving 
contributions of order p", and the sum has an upper cutoff 
at n = j, since ry corresponds to the jth generation singlets. 
Assuming that a„ = a\U"~ l , for some constant a (which can 
be numerically checked to be a reasonable approximation for 
small p), we have 



1 +fliP + fl2P 



i , 1 - («P) ; ' 

-ajp j = 1 + — -flip, 

1 - ocp 



and taking into account that ry ~ roZ J sm we can write 

(ap)'~r^><» with T1 (p) = i-^P). 
Combining the above results we conclude that 

C*(r,)~l 

for p < oc~ 1 , reproducing the zeroth order MDH prediction, 
but a nonuniversal behavior 



C xx (rj) 



1 



Jl(p) 



is obtained for p > oc~ 1 . 

For the silver-mean XX chains, an estimate of the a„ based 
on numerical implementations of the MDH method gives 
a i ~ 9.5 and a ~ 3, but with some dependence on ry and p. 
As shown in Fig. |SJ the decay exponent r| xx of the average 
correlations approaches unity as p — > 0, but starts to decrease 
more rapidly for p > 0. 1 , considerably less than 1 /a; this dis- 
crepancy indicates that Eq. d!6l >. with the assumption of con- 
stant a„ 's, although providing a valuable insight into the origin 
of the nonuniversal behavior, is not a good approximation for 
larger coupling ratios. The exponents predicted by the second- 
order MDH scheme are systematically smaller than the pre- 
sumably exact ones obtained from the free-fermion method 
(which should tend to r\ xx = 1/2 as p — > 1). This also happens 
for the decay exponent r\f x p of the typical correlations, which 
diverges as p — > 0, in agreement with the fact that, in this limit, 
the chain decomposes into independent singlets. A similar be- 
havior is observed for the transverse correlations C K (r) and 



Tltyp = 2 



C t Z y p (r) (but now the decay exponents approach r|" 
as p — > 1). 

On the other hand, dominant ground-state correlations in 
the Heisenberg silver-mean chain closely follow the predic- 
tions of the zero th-order MDH scheme, as can be seen in Fig. 
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Figure 8: Decay exponents of C xx (rj) and Q" p (r) for the XX silver- 
mean chain as a function of the coupling ratio p = J a /Jb, obtained 
from both the second-order MDH scheme (chains with 8819 sites) 
and the free-fermion method (chains with 3363 sites). Errors bars 
are at most the same size as the symbols themselves. 



0b). This is due to the fact that the effective coupling ratio 
decreases as the RG proceeds, and the contribution to C(r ; ) 
due to spin pairs other than those connected by strong bonds 
becomes exponentially negligible. Typical correlations decay 
not as a power law, but rather according to 



j>o 0--C O-' 



Ctyp(r) 



-«ypln 2 (r/r ) 



precisely the same form of the dynamical scaling; by fitting 
the numerical results, the constant /j typ is found to be approx- 
imately fi/2, with fj given by Eq. il5\ . As in the random- 
singlet phase, the scaling form of the typical correlations is 
similar to that of the lowest gaps, reflecting the fact that two 
spins separated by a distance r are basically uncorrected until 
the energy scale is of order A(r), when they become weakly 
correlated through an intervening spin taking part in a singlet 
pair. 



C. The bronze-mean sequence 

The bronze-mean sequence is built from the substitution 
rule a — > aaab, b — > a, with a large XX rescaling factor 17 



'3 + VT3" 



36.03. 



The bond-distribution attractor produced by the MDH RG 
scheme is not a fixed point, but a two-cycle; 64 apart from a 
few bonds near the chain ends, the same distributions alter- 
nate between even and odd generations, as shown for J a < Jb 
in Fig. [9] In this case, the second-generation couplings relate 
to the original couplings by 



J 

h = *kh-k and & = 



' / 4 
J b 



4' 



Figure 9: First five generations of the bronze-mean XXZ chain with 
Ja < Jb, eacn showing the leftmost 24 sites. The numbers indicate 
the positions of the sites in the original chain. Encircled blocks con- 
tribute effective spins when renormalized. The labels on the right 
denote the effective bonds in each generation. The attractor of the 
block distribution is a two-cycle, reached at the second generation. 
All bonds to the right of the horizontal arrow follow the same se- 
quence in the second (third) and fourth (fifth) generations. 



Likewise, in terms of the couplings in the previous generation, 
we write the third-generation couplings, 

p 

/4=73Y4-r and J B = %J a , 

Jb 

and the fourth-generation couplings, 

/>7t7f| and f h 

Since the attractor is now a two-cycle, and not a fixed 
point, we must relate the couplings in the fourth and second- 
generations. By eliminating Ja and Jg in the above equations, 
we get 

J b J b 
so that the coupling ratios satisfy the recursion relation 

/' / 



Jb 



while the corresponding gaps are related by 

A' V 



Jb 



: Y2Y3^P 4 



The distance between spins connected by strong bonds in the 
jth generation is rj = 1, 13, 43, 469, 1549,. . ., which asymp- 
totically gives rj+2/fj — %m> so that r^j ~ roi ] hm . Thus, solv- 
ing the recursion relations for p and A, and taking into account 
that p = 72 p, we obtain the dynamic scaling form 



^—^p'expf-^ln 2 ^ 



(17) 
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Figure 10: Ground-state correlation functions of the XX (a) and 
Heisenberg (b) bronze-mean chains, obtained from the zeroth-order 
MDH scheme (chains with 2010601 sites). 



with 



C(P) 



in (-yfyfytP 4 

lnTbm 



and 



In (Y2Y4 
ln 2 T bm 



Of course, the same form is obtained if we choose to look 
at the odd generations. In the XX limit, as 72 = 74 = 1, we 
again have /j = and Eq. (I17> corresponds to a nonuniver- 
sal power-law scaling behavior, with a dynamical exponent 
z = C(p); once more, as in all marginal XX chains, £ equals 
the leading term in the p <C 1 expansion of the exact dynam- 
ical exponent^ As in the Fibonacci and silver-mean chains, 
choosing J a > //, leads to the same scaling behavior, since af- 
ter the first lattice sweep the bond distribution is essentially 
equal to the one obtained for J a < Jy. 

Thus, the bronze-mean chains present qualitatively the 
same low-temperature thermodynamic behavior as the Fi- 
bonacci chain. However, this is not the case for ground-state 
properties. As indicated in Fig. [5] the renormalization process 
in each generation involves all spins in the chain, and gives 
rise to a hierarchy of effective spins, analogous to that shown 
in Fig. ^2 As a consequence, an effective spin in the y'th 
generation represents 3 7_1 real spins. So, while the ground 
states of the Fibonacci and silver-mean chains could be de- 
scribed as 'aperiodic singlet phases', from which excitations 
of a given energy involve spins separated by a single, well- 
defined distance, low-energy excitations in the bronze-mean 
chain involve an exponentially large number of spin pairs, 
whose distances are distributed in an increasing range. This is 
reflected in the ground-state correlation functions, which ex- 
hibit a fractal-like structure, as seen in Fig. ^| The strongest 
correlations in the chains correspond to the distances rij = 13, 
469, 16897, . . ., and their scaling behavior can be obtained by 
the following analysis. 

Consider a pair of neighboring effective spins belonging to 
the same block in the /'th generation, and let cq be their zeroth- 
order correlation. Each of these spins represents 3-*~ l real 
spins, so that for each such pair there are 3 J ~' pairs of real 
spins separated by the same distance ri contributing to the to- 
tal correlation per site C(rj). However, the contribution of a 



real pair to C(rj) depends on the string of Clebsch-Gordan co- 
efficients indicating the weight of its two spins in the effective 
spins: each time the intermediate effective spin representing 
a real spin Sk is located at the ends (the center) of a three- 
spin block, the weight of Sk is multiplied by a factor C13 (02,3) 
upon renormalization. (These coefficients are in general dif- 
ferent for xx and zz correlations; see AppendixlAl) Since each 
effective spin in the jth generation has gone through j — 1 
renormalizations, a real pair can be classified according to the 
number n of factors ci 3 present in the (equal) weights of its 
spins. The contribution of all type-n pairs to C(r ; ) is propor- 
tional to the number of such pairs, being given by 



(") _ 2«_ 



2 \« , l \ ./-!-« 

" n\{j-\-n)\ K ~ lx 



! )"(4 3 ) J 



Thus, the total contribution of a single effective-spin pair to 

C{rj)h 



Sj 



L^ ) = (43+2 C ,, 3 ) ;W koh 



which gives 



C{rj)~n jg j 



where nj is the fraction of active spins in the jth generation. 
Since asymptotically we have j = ln(r ; /ro) / In -y/Tbm, this last 
result can be written as 



C(rj) 



with T| = 1 — 



ln(43 + 2c u) 



(18) 



For the Heisenberg chain, ci, 3 = 73 = 2 /3 and C23 = : /3, so that 
T| = 1 . For the XX chain, r| depends on whether we look at 
longitudinal or transverse correlations: in the former case we 
have cf 3 = 1/V2 and cf* 3 = 1/2, so that rf* ~ 0.875, while in 
the latter case cf 3 — 1/2, cf 3 = 0, and so T| zz ~ 1.387. These 
values are fully compatible with the results from numerical 
implementations of the MDH scheme shown in Fig. ^3 an d 
agree very well with free-fermion calculations for XX chains 
with p<l. Again, larger coupling ratios lead to nonuniversal 
decay of the correlations, except in the Heisenberg limit. 



D. A sequence producing effective-spin triples 

The appearance of an effective-spin hierarchy is better il- 
lustrated by the sequence obtained from the substitution rule 
a — > bbaba, b — > bba. The first three generations of the chains, 
for J a < Jt,, are shown in Fig. ^2 f° r Ja > Jb the first lattice 
sweep interchanges the roles of weak and strong bonds, recov- 
ering the former case. Renormalization involves both three- 
spin and two-spin blocks, and each lattice sweep reproduces 
the original sequence, yielding effective couplings given by 

J' a = Y2Y3 f and 4=T&. 
Jb 
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Figure 11: First three generations of the XXZ chain with cou- 
plings following the sequence in Sec. IV Dl for J a < showing 
an effective-spin hierarchy. Solid lines indicate strong bonds; for 
clarity, weak bonds are not represented. Shaded blocks contribute 
effective spins when renormalized, while white blocks form singlets. 
A third-generation effective spin represents three second-generation 
and nine first-generation spins. 
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Figure 12: First two generations of the XXZ chain with couplings fol- 
lowing the sequence in Sec. IV HI for J a < Jb- The numbers indicate 
the position of the spins in the original chain. The first lattice sweep 
generates a fixed-point bond distribution with four different effective 
couplings J a through Jj. 



so that the recursion relations for the coupling ratio and the 
energy gap are 

p' = 72 p and A' = y^pA. 

The size of three-spin blocks follows r, = 2, 6, 22,. . ., while 
that of two-spin blocks corresponds to rj/2, leading asymp- 
totically to a rescaling factor 



tst 



2 + V3-3.73, 



and a dynamical scaling relation 



with 



C(p) = - : 



In 



(y 2 1/2 YsP 
lnT st 



ln(yf 

and n = \ — 

lrrTst 



Again, aperiodicity induces nonuniversal behavior for < 
A < 1, and a weakly exponential scaling in the Heisenberg 
limit. 

As in the bronze-mean chains, discussed in the previous 
subsection, excitations of a given energy involve an exponen- 
tially large number of spins, due to the effective-spin hierar- 
chy. More precisely, since each effective spin in jth gener- 
ation represents 3 7_1 real spins, excitations with energy A,-, 
corresponding to breaking a y'th generation singlet, involve 
2 • 3 ; ~' spins; exciting a three-spin block in the same genera- 
tion costs an energy of the same order, and involves 3 ' spins. 
Dominant ground-state correlation functions also decay as in 
Eq. Jl 81 . but now with 



In 



11 = 1 



(43 



-2c? 3 



hlTst 



yielding rf* ~ 0.830 and r\ zz ~ 1.526 for the XX chain, and 
T| = 1 for the Heisenberg chain. These values are again fully 
compatible with results from numerical implementations of 
the MDH scheme. 



E. A marginal tripling sequence 

This sequence is generated by the substitution rule a — > aba, 
b — > bba. As discussed in Ref. 17, this type of aperiodicity 



may lead to marginal behavior even in anisotropic XY chains. 
As shown in Fig. [^] for J a < Jb the MDH scheme produces 
a second-generation lattice with four different effective cou- 
plings, given by 



J, 



^=Y2Y3T-, h = Y3 -fa, 
Jb 



Jc = %Jai h 



with an effective coupling ratio p = J a /Jb — J2P- (Choosing 
/„ > Jb interchanges the roles of f c and J c [, otherwise produc- 
ing the same bond distribution.) The bond distribution does 
not change upon further lattice sweeps, and the effective cou- 
plings satisfy the recursion relations 

Jb Jb 



7t Ja^c J a^d 

Jb Jb 

Noting that J a = JJd /Jb , we can write 

V 1 

J'b h 

so that aperiodicity is marginal even in the Heisenberg limit. 
The recursion relation for the gaps is 

A' = ^A = y 2 pA, 
Jb 

and the size of the singlets formed along the generations fol- 
lows rj = 3' — 1, with a rescaling factor T mt = 3, so that the 
dynamic scaling relation is given by 

j\j - rj 

with a nonuniversal dynamical exponent 

In (y2p) 



Z = C(P) = - 
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Figure 14: First two generations of the binary Rudin-Shapiro XXZ 
chain, for J a < J b . The first lattice sweep generates 8 different ef- 
fective couplings, J a through J/,, labeled in the figure by the letters 
a — h. Further renormalization does not change the bond distribu- 
tion. Starting from the second generation, only blocks coupled by J b 
(thick lines) and J c (thin lines) bonds are renormalized. (For clarity, 
weaker bonds are not drawn in the picture.) 



Figure 13: (Color online) Log-log plots of the specific heat c(T) and 
susceptibility %{T) as functions of temperature fox XXZ chains with 
couplings following the marginal tripling sequence, for J a /Jb = V 10 
and three different values of the uniform anisotropy A. 

Thermodynamic properties can be estimated by using the 
same idea of the 'independent-singlet' approximation de- 
scribed for Fibonacci chains, with slight modifications due 
to the fact that the first lattice sweep (but not the later ones) 
involves renormalization of both two- and three-spin blocks. 
Thus, the free energy per site can be calculated by adding to 
Eq. d 1 3I > a term representing the contribution of spins renor- 
malized in the first lattice sweep, and given by 

fi{h,T) = l -F vm ,{J b ,A h \hJ)+ l -F t ^ Xe {J b ,A h \hJ), (19) 

where ^triple {J,A;h,T) is the free energy of a spin triple obey- 
ing the Hamiltonian 

^triple = /(SjSI + SjSj + ASjS^) 

+ /(SfSf +S y 2 Sl+AS z 2 S z 3 ) 
- h(S\+S z 2 + Sl). 

Notice that three-spin blocks yield effective spins when renor- 
malized, and these will pair with other real or effective spins 
to form singlets in the second lattice sweep, but this is not 
taken into account by Eq. dl9l >. In order to obtain a correct 
estimate of the low-temperature susceptibility, we must mul- 
tiply the contribution arising from spin triples by a factor like 
e~ Jb ' T . Results for the temperature dependence of the spe- 
cific heat and susceptibility are shown in Fig. [D] for p = i/io 
and three values of the uniform anisotropy A, corresponding 
to the XX and Heisenberg chains and to an intermediate case. 
Both quantities exhibit log-periodic oscillations, obeying the 
scaling forms 

«m - zcn~^'<%(|£). 

with A. being the asymptotic ratio between the gaps in suc- 
cessive generations, while G c and G z are periodic functions 
(with period one). In the XX and Heisenberg limits, we 
have A, = y^p; for intermediate anisotropies, A equals p e g-, a 
coupling ratio defined at the energy scale in which effective 
anisotropies become negligible. 
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Figure 15: Pattern of three-spin blocks and isolated spins leading 
to the effective-spin hierarchy in Rudin-Shapiro XXZ chains. Thick 
lines indicate strong bonds. Shaded blocks contribute effective spins 
when renormalized; white blocks form singlets. 



Also as a consequence of the strictly marginal character of 
aperiodic fluctuations for all anisotropies in the regime < 
A < 1, dominant ground-state correlations follow C(r,-) ~ 1/r; 
in the p -C 1 regime, but nonuniversal behavior should be ob- 
served for larger coupling ratios. 



VI. RELEVANT APERIODICITY 

A. The binary Rudin-Shapiro sequence 

The Rudin-Shapiro sequence is originally defined as a four- 
letter sequence, 53 generated by the substitution rule a — > ac, 
b — > dc, c — > ab and d — > db. It has the interesting property 
that its geometrical fluctuations mimic those induced by a ran- 
dom distribution. In order to reduce it to a binary sequence, 
we make the associations c = a and d = b, obtaining an in- 
flation rule for letter pairs, given by aa — > aaab, ab — ► aaba, 
ba — > bbab and bb — > bbba. The rule generates blocks having 
between 2 and 5 spins, and is symmetric under the interchange 
of a and b, so that the scaling behavior is invariant with respect 
to the interchange of J a and J b . The left-end of the first two 
generations of the Rudin-Shapiro chains is shown in Fig. \^\ 
for J a < J b . 

Blocks with more than 3 spins are eliminated in the first lat- 
tice sweep and do not appear in later generations. Both two- 
and three-spin blocks are present in the fixed-point block dis- 
tribution (already reached at the second generation), and upon 
renormalization the sequence produces an effective-spin hier- 
archy, stemming from approximately mirror-symmetric pat- 
terns of three-spin and five-spin blocks in the original lattice. 
This is illustrated in Figs. H5landfT6l 

In the jth generation, three-spin blocks have size rj = 
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Figure 16: Pattern of five-spin and three-spin blocks leading to the 
effective-spin hierarchy in Rudin-Shapiro XXZ chains. Thick lines 
indicate strong bonds. Shaded blocks contribute effective spins when 
renormalized; white blocks form singlets. 



2 • 4 ; ~ 1 (with a rescaling factor T ls = 4), while two-spin blocks 
have size rj/2. The first lattice sweep generates effective cou- 
plings J, having 8 different values, 




/«=Y2Y3Y44, fb = 1375 Ja, 



J i i 

Jc = J2J3 -f ' & = WW T" ' J' = "fc™ ~f~ ' 
Jb Jb Jb 



Ji 



Ji 



/,-*»4 /.-**,& 



Ji 



r 



and whose bond distribution remains unchanged upon renor- 
malization, leading to the recursion relations 



p 3^aJ~dJg ,-/ 

J a —ii—rT, 



J~bJ} 



J b =J3Jf, 



Jb Jc 



ft .2., JaJdJf p 2 JqJdJg p ^JeJfjf, 

Jf = Y2Y3 f n ' J g Y2Y3 r f , Jh = Y2- 



JJi 



JbJc 



Pi 

J b Jc 



Defining a new effective coupling Jo — JaJd/Jo we obtain a 
system of three recursion relations, 



f _ v 4 JoJf 
J — Y2Y3 — p~ ; 

J b 



f b =f3ff, and T f = $t 



JoJf 
Jb 



With coupling ratios p = Jo/Jb and a = Jf/Jb, and a gap pro- 
portional to ft,, we have 



Pj=l2p 2 j-U Gj=y2P~j-l, and —J-=y 3 Cj, 



which after eliminating 6j yields 



Figure 17: Ground-state correlations for chains with Rudin-Shapiro 
couplings, obtained from extrapolation of numerical MDH results for 
chains with 2 16 to 2 20 sites, (a) C xx {r) (upper solid curve) and C zz (r) 
(lower solid curve) for the XX chain. The dotted curve is proportional 
tor -3 / 2 . (Curves offset for clarity.) Inset: dominant C" correlations, 
corresponding to distances rj = 2- 4' , fitted by a law of the form 
rC xx (r) = yo+yi lnr (dashed curve), (b) C(r) for the Heisenberg 
chain (solid curve). The dotted curve is proportional to 1 jr. 



with CO = 1/2, 



c= ln(WY|) and 



lnx r . 



-ln^p 2 



So we obtain, for the whole regime < A < 1, the dynamical 
scaling form predicted for the XX chain, reproducing the result 
for the random-singlet phase. 

For chains with RS couplings, effective-spin formation de- 
termines the dominant ground-state correlations, but the cor- 
responding hierarchy is slightly different from the ones seen in 
sees. IVClandlVTJl now involving both three-spin (and some 
five-spin) blocks and unrenormalized spins. As illustrated in 
Figs. EJand^J for each block renormalized in the j'th gener- 
ation the correlation between its end spins connects a number 
of order V original spin pairs separated by the same distance 
rj (the size of the block), yielding a contribution to the aver- 
age correlation in the Heisenberg chain and C xx (rj) in the XX 
chain given by 



8i 



2c 



1,3 J 



i 3 1 - 

k=l 



2c 



1.3 



ko|, 



where co is the correlation between end spins in a three-spin 
block. For the Heisenberg chain 2c\ 3 = 8 /9 < 1, and thus 



C(rj)~ gj . nj , 



1 



(20) 



P>=Y2P;-i and 



A 7 ~ l /2 

A/-1 



Solving the recursion relations we obtain 



r 

Aj ~ rj exp 



I- 1 1 - 

JO 



where rij ~ 1/ rj is the fraction of three-spin blocks in the jth 

generation. For the XX chain 2 fc^O = 1, so that gj has a 

term proportional to j, and C xx (rj) carries a logarithmic cor- 
rection, 



(21) 



15 



Figure 18: First two generations of the 6-3 sequence, discussed in 
Sec. IV1BI for J a < J{,. In the second generation, dashed lines in- 
dicate effective J a bonds, while thick and thin solid lines denote ef- 
fective J/, and J c bonds. In subsequent generations, the J c J a pairs 
change to J a J c , and the first J a becomes J c , but the bond distribution 
is otherwise unchanged. 



where yo and y\ are constants. The zz correlation between 
end spins in a three-spin block is zero, so that the dominant 
correlations correspond to spin pairs (connected through one 
of the effective end spins and the middle spin) at distances 

= 1/2. 



r) = ±\}-t ±\i-i ± • • • ± 1, with average 



whose contribution is given by g'j 
We then have 



f-l 



since c 



1.3 



-3/2 



(22) 



Eqs. \2\\ and J22l i should be contrasted with the random- 
singlet isotropic result C{r) ~ r~ 2 , indicating a clear distinc- 
tion between the ground-state phases induced by disorder and 
aperiodicity, even in the presence of similar geometric fluctu- 
ations. This is related to the inflation symmetry of the aperi- 
odic sequences, which is absent in the random-bond case (or 
in aperiodic systems with random perturbations^ 4 -). Its effects 
are exemplified by the fractal structure of the ground-state 
correlations visible in Fig. which displays results from 
numerical implementations of the MDH method for both XX 
and Heisenberg chains, showing conformance to the scaling 
forms in Eqs. J20l- (l22i . Contrary to the marginal sequences, 
these scaling forms should be observed in the large-distance 
behavior of Rudin-Shapiro XXZ chains for any coupling ratio 
p 7^ 1 ; we expect a crossover from the uniform to aperiodic 
scaling behavior as larger distances are probed for p close to 
unity. Free-fermion calculations in the XX limit support this 
picture. 



Thus, defining the effective coupling ratios 



_ J a 7 | J a 



i ~ Jc Ja 

and a = — =j2 — , 
h h 



we can rewrite the recursion relations as 

A' /' 

P'=^P 2 , 6' = Y2p, and — = ^=y 2 dp. 
In the jth generation, we have a 2 = Pj, and thus 

P;+i=Y2Pi and A ;+i =Y2pJ 2 A ; - 

The length of the singlets correspond to rj = 1, 9, 45, 225,. . ., 
so that asymptotically rj ~ ryX-', with ro = 9 /25 and x = X63 = 
5. Solving the above recursion relations we obtain the dynam- 
ical scaling behavior, 



r j exp 



(23) 



with a wandering exponent co = In 2/ In 5 ~ 0.431 and 



31ny2 
In 5 



and 



A/=-^ln(Y^p 3 ) 



where p = J a /Jb is the original coupling ratio. 

If we choose J a > blocks with 2, 3 and 4 spins cou- 
pled by strong bonds appear along the chain. Effective spins 
are produced by the first lattice sweep, yielding effective cou- 
plings 



Ja 



j3 

Jb= Y2Y3Y'' b 



7 2 ' 



and / C =Y2^^> 



whose distribution is the same as that of the third-generation 
bonds for J a < J/,, and which remains unchanged upon renor- 
malization. Thus, the scaling behavior is the same as above, 
but now with a 'bare' coupling ratio 



B. The 6-3 sequence 



Ja 7 7 ( Jb 



This sequence is generated by the substitution a — > babaaa, 
b — > baa, and its XX wandering exponent is CO = In2/ln5, with 
a rescaling factor X63 = 5. Application of the MDH scheme 
leads to a fixed-point bond distribution with singlet renormal- 
ization only, so that no effective-spin hierarchy is present. For 
J a < Jb, as depicted in Fig. ^] three effective couplings are 
produced after the first lattice sweep, 



Ja 



-T2 



■11 

I 2 ' 

J b 



Jb =Ja; Jc = j2 



J l 

Jb' 



and upon further lattice sweeps we obtain the recursion rela- 
tions 



Thermodynamic properties can be estimated as in the Fi- 
bonacci case, by using the 'independent-singlet' approxima- 
tion. Plots of the specific heat c(T) and susceptibility %(T) 
as functions of temperature are shown in Fig. ^] and com- 
pare quite well with results from numerical diagonalization, 
even for relatively large coupling ratios {J c ,/Jb = V 4 )- This is 
not surprising, given the fact that the effective coupling ratio 
rapidly decreases as the RG proceeds, even for the XX chain. 
As seen in the inset of Fig. I19t b). at temperatures of the order 
of the gaps A ; the susceptibility follows the scaling form 



1 



rllnT 



l/co' 



Pi 



■Ti—fi 



Jl 



J Ore 

h 



pi 

J a J c 
J b 



which can be readily obtained from Eq. d23t by assuming 
that singlet pairs are magnetically frozen, while active spins 
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Figure 19: Thermal dependence of the specific heat (a) and suscepti- 
bility (b) of the XX chain with couplings following the 6-3 sequence 
of Sec. I VI Bl obtained for J a /Jb = V 4 fr° m both numerical diago- 
nalization of chains with 46875 sites (circles), and the MDH scheme 
(solid curves). The inset in (b) presents a log-linear plot of T versus 
(xT)~ (i> , with CO = In2/ln5, showing that for T ~ A;, corresponding 
to the specific-heat maxima, the susceptibility satisfies the scaling 



formx~r _i |lnr 



-1/to 



(dashed line); at intermediate temperatures, 



a Curie-like behavior is observed. 




Figure 20: (Color online) Ground-state correlations of the XX chain 
with couplings following the 6-3 sequence of Sec. lVlBl for two dif- 
ferent values of the coupling ratio p = J a /Jb, as obtained by numer- 
ical diagonalization of chains with 1874 sites. Peaks in the curves 
correspond to the characteristic distances rj = 9, 45 and 225. 



contribute Curie terms to %. Estimates of c(T) and %(T) for 
chains with anisotropies < A < 1 are qualitatively identical 
to the ones for the XX chain. 

Since no effective-spin hierarchy is present, and aperiod- 
icity is relevant, dominant ground-state correlations, for any 
coupling ratio p ^ 1 and sufficiently large characteristic dis- 
tances rj, should decay as 

c«( )~c-(o)~-, 

for all anisotropies in the regime < A < 1 . This is confirmed 
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Figure 2 1 : Leftmost 24 sites in the first six generations of XXZ chains 
with couplings following the fivefold- symmetry sequence, discussed 
in Sec. IVICI for J a < J/,. The attractor of the bond distribution is 
a two-cycle, reached after three lattice sweeps. Second-generation 
bonds are denoted by J a through /g; third- and fifth-generation cou- 
plings are labeled a through /, while A through G label fourth- and 
sixth-generation bonds. Couplings Jq only occur much farther along 
the chains. Starting from the fourth generation, lines indicate blocks 
to be renormalized. 



in the XX limit by numerical diagonalization, as shown in Fig. 



C. The fivefold-symmetry sequence 

The sequence produced by the substitution rule a — * aaab, 
b — > bba is related to binary tilings of the plane with fivefold 
symmetry. 36 The (quite large) XX rescaling factor is % = 25 + 
10\/5 ~ 47.36, with a wandering exponent co = In 3/ In ~ 
0.285. 

Under a numerical implementation of the MDH scheme 
with J a < Jb, we obtain a quite intricate pattern: after a four- 
bond transient produced by the first two lattice sweeps, a two- 
cycle periodic attractor is reached, where six- and seven-bond 
distributions alternate, as depicted in Fig. |^ (With J a > Jb, 
the same two-cycle is reached after the first lattice sweep.) 
The distance between spins connected by the strongest bonds 
in each generation correspond to r\ = 1, 3, 33, 190, 1385, 
9050, which asymptotically gives r/+2/V; ~ is ■ The 
equations relating the effective couplings of the fourth and 
fifth generations are 



Ja 



3 J\JcJ~F 



J~b = j3J~E, Jc = yi 



JaJcJg 

p ' 
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Jf = 727: 
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while between the couplings of the fifth and fourth generations 
we have 

/1=7^7 3 #, 4 = ™^£, fc = ^, 

JbJc Jb Jb 
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D — 



2 ^3 fTf i J E 



j' f =U 
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Eliminating the fifth-generation couplings and defining the 
ratios 



- J~A 

9 = T B 



01 = v 



2 



P 
J F 

JaJb 



and 03 = — , 
Je 



we can write a set of four recursion relations, 

P' =72P>, 



73 



o 2 = 72 Y3P (01O3) , 03=75730103. 
The gaps in successive even generations obey 

A ' J'n s 2 
X = ^=T 2 73Pa 1 a 2 , 

and by expressing 01,02 and 03 in terms of p we get 

where now 2j + 2 labels the lattice generation and a is a con- 
stant depending on the values of the coupling ratios in the 
fourth generation. Solving this last equation gives 



,\7>" C 



3.' 



with 



A=ar 2 \ B=tf\ and C = Jf~9? 



For large enough /', since j = ln(r//Vo)/lnTff, we have 

(O-l 



exp 



with 



/j = — InC and co = 



In 3 
lnTff ' 



again obtaining, for the whole anisotropy regime < A < 1, 
the same scaling form predicted for the XX chain. 

The effective-spin hierarchy produced by the RG process is 
analogous to that in the bronze-mean chains, so that ground- 
state correlations behave as in Eq. |T8) , with r| = 1 in the 
Heisenberg chain, but Tf* ~ 0.884 and T[ a ~ 1 .359 in the XX 
limit. However, these figures are not so well reproduced in 
the numerical calculations, even for chains with ~ 1 .6 x 
10 6 sites, most probably due to the extremely large rescaling 
factor. 



VII. DISCUSSION AND CONCLUSIONS 

For all aperiodic sequences discussed in the previous sec- 
tions, the recursion relations for the main coupling ratio and 
the energy gaps have the forms 



Pi+i= c Pj and A i+i = /1/2 PyA/' 



(24) 



where c, f\ and fa are A-dependent nonuniversal constants, 
and I (a rational number) and k (an integer) relate to the 
number of singlets involved in determining the effective cou- 
plings. In particular, k is ultimately the difference in the num- 
ber of singlets producing the effective couplings whose ratio 
is p < 1. 

If k > 2, the recursion relation for p always has a stable 
fixed point at p* = 0, so that the effective coupling ratios 
become exponentially small as the renormalization proceeds, 
indicating that asymptotic results obtained from the MDH 
method should be essentially exact. Taking into account the 
scaling behavior of the characteristic distances rj ~ roT 7 , Eqs. 
(1241 lead to the dynamical scaling form 



-C Vln 2 (r ; /r ) - M (rj/r )<» _ -p( rj /r f 



(25) 



with £ and // nonuniversal constants, 



In 



k(k-l) 



p being the 'bare' coupling ratio, and 



co = 



Ink 
InT 



Note that co has the same form as the exact wandering expo- 
nent for XX chains with nondimerizing aperiodic couplings, 
given in Eq. Q. Moreover, CO depends only on the topology 
and the self-similar properties of the sequence, being indepen- 
dent of the anisotropy in the regime < A < 1 . 

If k = 1, the recursion relation for p has a line of fixed 
points, provided that c = 1, which is generically the case in 
the XX limit; otherwise p* = is a stable fixed point. The 
general solution to Eqs. i24\ is 



(26) 



where 



C(p) = - 



ln(/./ 2 'VV 



lnx 



and fi = — - 



In 



,Vy/ 2 



ln 2 T 



Unless fa^\, which, among the sequences studied here, hap- 
pens only for the relevant fivefold-symmetry sequence of Sec. 
IVICI /j is zero in the XX limit. This means that we can identify 
a nonuniversal dynamical exponent z = £(P)> and the scaling 
behavior of thermodynamic properties depends on the cou- 
pling ratio for the whole anisotropy regime < A < 1. In 
the Heisenberg limit (A = 1), unless /j — 0, as in the marginal 
tripling sequence of Sec. IV El Eq. i26\ describes a weakly 
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exponential dynamic scaling. In this case, aperiodicity can 
be viewed as a marginally relevant operator (co — * + ) in the 
renormalization-group sense. 

These results strongly suggest that low-temperature ther- 
modynamic properties of any antiferromagnetic XXZ chain 
with anisotropies intermediate between the XX and Heisen- 
berg limits, and couplings following a given binary aperi- 
odic sequence, can be classified according to a single wan- 
dering exponent CO, which is known exactly for XX chains. 
This generalizes what happens in random-bond XXZ chains 
(for which CO = 1/2), where thermodynamic properties in the 
anisotropy regime — !/2 < A < 1 are those characterizing the 
random-singlet phase. 8 - 41 Note that, although the above clas- 
sification seems to imply an anisotropy-independent critical 
value co c = for the relevance of aperiodic fluctuations on 
the low-temperature behavior of XXZ chains, it does not show 
that CO plays the role of a genuine wandering exponent, in the 
sense that fluctuations scale as g ~ N a , for general easy-plane 
anisotropies. In any case, due to the fact that the critical expo- 
nents (including the correlation-length exponent v) of the uni- 
form XXZ chain are known to vary with the anisotropy along 
the whole critical line — 1 < A < l^S^L it remains an open 
question how the present results fit into the framework of the 
Harris-Luck criterion. 

Of course, Eqs. (I24> are valid for all anisotropies < A < 1 
only if the bond distribution generated by the MDH method is 
independent of A. This is certainly the case for strong enough 
modulation. (How strong this modulation has to be depends 
on the various block sizes produced by the sequence.) How- 
ever, from numerical implementations of the method, we find 
that, even when the blocks selected for renormalization in the 
first few lattice sweeps depend on A, a universal distribution is 
eventually reached, in much the same way as when we choose 
J a > Ji, instead of J a < Jt,. Thus, we expect that, for general bi- 
nary substitution rules inducing relevant aperiodicity, the scal- 
ing form in Eq. d25t holds for all coupling ratios p 7^ 1 . 

An approximate picture of the ground state and of the low- 
est excitations in the presence of aperiodic couplings can 
also be deduced from the MDH scheme, and is revealed by 
the behavior of the pair correlation functions. As the en- 
ergy scale is reduced, two types of behavior can be iden- 
tified: either the RG process produces a hierarchy of sin- 
glets (as in the Fibonacci, silver-mean, marginal-tripling, 
and 6-3 sequences), or a hierarchy of effective spins (as in 
the bronze-mean, spin-triple, Rudin-Shapiro, and fivefold- 
symmetry sequences). The first type reveals a kind of self- 
similar, 'aperiodic-singlet' phase, from which (singlet-triplet) 
excitations involve strongly coupled pairs composed of spins 
separated by well-defined characteristic distances. In the sec- 
ond type, since the number of spins contributing to an effec- 
tive spin increases exponentially along the hierarchy, excita- 
tions of a certain energy involve spins separated by a wide 
range of distances, giving rise to a fractal structure of the cor- 
relation functions. Notice that, contrary to the finite temper- 
ature behavior, there is no relation between the ground-state 
properties and the marginal or relevant character of the aperi- 
odicity. 

For aperiodic sequences inducing strictly marginal fluc- 



tuations, we could account for the nonuniversality of the 
correlation-function decay exponents by a numerical calcu- 
lation based on a second-order expansion of the ground-state 
vectors. This compares quite well with results from numerical 
diagonalization in the XX limit, which indeed show that the 
zeroth-order MDH predictions are reproduced in the strong- 
modulation regime. 

The results on relevant aperiodic couplings show that ge- 
ometrical fluctuations, measured by the wandering exponent 
CO, are not determinant for ground-state properties, although 
they control the low-energy scaling of thermodynamic quan- 
tities. In particular, both random bonds and Rudin-Shapiro 
couplings are characterized by co = 1/2; however, correlations 
in the random-singlet phase are entirely different from those 
in XXZ chains with Rudin-Shapiro couplings. This is a con- 
sequence of the inflation symmetry induced by substitution 
rules, which is clearly absent in random chains. (Analogously, 
comparative studies 55,56 between random-bond and Rudin- 
Shapiro quantum Ising chains show that, although the corre- 
sponding scaling properties are similar at the critical point, 
only randomness is capable of producing the off-critical Grif- 
fiths singularities!^ 2 * 5 .!) Nevertheless, aperiodic and random 
XXZ chains share the feature that average and typical behav- 
ior are strikingly distinct, and that average correlations de- 
cay as power laws. Finally, aperiodic ground-state phases 
are unstable towards random perturbations, which break infla- 
tion symmetry, and the random-singlet behavior is ultimately 
recovered^ 
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Appendix A: RENORMALIZATION OF MULTI-SPIN 
BLOCKS 

In this Appendix, we derive the expressions for the renor- 
malized coupling constants used in the extension of the Ma- 
Dasgupta-Hu method to XXZ chains with aperiodic couplings. 

Contrary to the random-bond chains discussed in Sec. |H] 
when couplings follow aperiodic sequences generated by in- 
flation rules we generally need to consider spin blocks with 
more than one strong bond, and thus more than two spins. 
For instance, in the Fibonacci sequence with J a > Jj, (see 
Fig. |2} there appear blocks with one or two J a bonds. Since 
we assume that all couplings are antiferromagnetic, the local 
ground state is a singlet for blocks with an even number of 
spins, but a doublet if the blocks contain an odd number of 

CO 

spins. 

Let us consider a block with n spins Si through S„ con- 
nected by equal bonds Jq, with anisotropy Aq. This is de- 
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Figure 22: Renormalization step involving the decimation of a two- 
spin block. 



scribed by the local Hamiltonian 



n-\ 



Hq=JoY, ( s r s i+i)Ao' 

where we introduced the notation 

(S ; -s ( i A s;s; ■ s]s) ■ .\s : s). 

The gap Ao between the ground-state energy of the block and 
its lowest excited multiplet depends on Jq and Aq. For two- 
spin and three-spin blocks we have 



A ( 2 ) = I+Aq^ ^ A ( 3 ) = 1 ( + ^ 8 j 



We define the strongest bonds in the chain as those producing 
spin blocks with the largest gaps Aq. 



An n-spin block to be renormalized is connected to its 
neighboring spins Si and S r by weaker bonds J\ and J r . The 
relevant part of the chain Hamiltonian is 



with 



H = H +Hi n 



H lr =J l (S r S 1 ) A[ +J r (S n -S r ) Ar 



(Al) 



The idea of the MDH method is to obtain recursion relations 
for the couplings by treating Hi, as a perturbation to Hq. 

If n is even (as in the two-spin case shown in Fig. 122b - 
the ground state of Hq is a singlet, which we denote by 
l^o), with a corresponding energy Eq. In the space of this 
singlet, the states of Si and S r are arbitrary. In the space 
spanned by the eigenstates |\P,-) of Hq (with energies E[) and 
the states \mi,m r ) = |m/) <g> \m r ) of St, r {mi r = ± l /i), the 
states \g(mi,m r )) = \mi,m r ) (g> |*Po) are degenerate. The first- 
order perturbative corrections to the ground-state energy Eq 
are zero, but the second-order corrections are given by the 
eigenvalues of the matrix 



V i i — 

r mj,m r ;m^,in r 



(g(mi,m r ) \Hi r \e) (e|fl/ r |g(wj,mj.)) 



Eo - Ei 

where the summation spans the excited states \e) = 
|m",m") ® (i= 1 , . . . , 2" — 1). In terms of the raising and 
lowering operators S^ 1 = S x ± iS y we have 



1 



>A=2\ S ' S J +S < S j) +AS i S P < A2 > 
and a little algebra shows that 



J 



V m/ ,m r ;m;X = J J l J r { m h m r\St S r \m' h m' r )Y, 



, , . ^ <^ \st | m*,) is+i * ) + opo \s+m |sr I *o> 



£ - £■/ 



1 / , c - c+ | , , v y (gojgfl ^) C% jgT I ^o) + QPo jg[ I gj) K | go) 
-J t J r (mi , m r | S, SJ | m ; , m r ) ^ 



Eq - Ej 



ami/ lc*e*l / 'W (^olSfl^X'J'/^l^ + ^ol^l^X^IS^o) 



¥0 



Eq - Ei 



r 



(A3) 



Since the first two terms on the right-hand side of Eq. ( IA3I are SU(2) symmetry is recovered). The above matrix elements 
complex conjugates, and noting that Eq — Ej is proportional to are precisely the ones corresponding to the Hamiltonian 
Jq, we can write 



Jll 



Vmw^m'r = Yn-^{mi,m r \S x l S x r + S y l Sl\m' l ,m' r ) 



JlJr 



+ y„ 8„ A/ A,- — - (mi,m r \ SjS z r \ m\ , m\. ) , 
where y„ and 8„ depend on Aq (with 8„ = 1 for Aq = 1, where 



H' = f (SJS X + Sffl + A'Spl) , 
with the effective parameters J' and A' given by 

/ = y n {ill an d a' = 8„A/A r . 
JQ 



(A4) 



20 



J,,A t Jo, A j n A r 

t / M \ 

5; S] S2 S3 S,- 




Figure 23: Renormalization step involving a three-spin block. 



For two-spin blocks we have 



Y2 



1 



1+A 



and 82 = 



1 



For larger blocks the parameters can be evaluated numerically 
as a function of Ao; however, for four-spin blocks we can ana- 
lytically determine 74 = 1 in the XX limit and 74 = 2/3 — V5/6 ~ 
0.378 in the Heisenberg chain. 

If n is odd (as in the three-spin case shown in Fig. 1231 . 
Hq has two degenerate ground states, which we denote by 
These can be associated with an effective spin- 1/2 
So, whose states can be described by the azimuthal quan- 
tum number mo, so that |mo = ±'/ 2 ) = |^o )• in tne s P ace U 
spanned by the states of So, S; and S r , the states |m; , m r , mo) 

(gi \m r ) ® |mo) are degenerate. The introduction of H\ r lifts 
this degeneracy, and we expect that, to order /; r /Jo< perturba- 
tion theory leads to an effective Hamiltonian H', with matrix 
elements given (apart from a constant) by 

H[ 111 =(mi,m r ,m()\Hi r \m', 1 m'm'n). 

Restricting ourselves to the space H, it is possible to write 

H = Jj (Si -So) A / +J r (So -S r ) A ^ , 

provided 

/ I (S f -Si) A( =^(S r So) A j (A5) 



and 



J r (S„ • S r ) Aj — J' r (So • S r ) A , . 



We now expand Eq. (IA5> with the help of Eq. (IA2> . and 
notice that 

(m/,m r ,m \SfS \ m,,m r ,m Q ) = 8,„/, m/ _ 1 8 m / imr 8 m / ) ^ +1 , 



(m,,m r ,m \S^S 1 | m' h m' r ,m' Q ) = S^^^jS^^ (m |Sj | m ) , 



(mi,m r ,m Q \S Z S^\ m',,m' r ,m' ) = m/m S m / ^S,,,/,,,,^ 



(mi,m r ,m \SfS\ \ m' h m' n m' ) = m/8 m ; ,,,,,8,,,; ^ (m |Sf |m ) , 

8, ; being the Kronecker symbol. By the Wigner-Eckart theo- 
rem, we can write 

(mo|Sf|m ) =Y„8 m j )imo+1 , 



<m |Sij I m ) = (8„y„) m 8 m /. mo , 

with y„ and 8„ depending on Ao, and we formally obtain the 
renormalized parameters 

J'l = Inh and A', = 8„A/ . (A6) 
Analogously, by symmetry we have 

j' r = y„J r and A;. = 8„A r . (A7) 
For three-spin blocks we obtain 

^Ao +v /a|78) 



Y3 = 



2+i (Ao + ^/Ag + sy 



and 



83 = t ( Ao + \/A 2 - 



while for larger blocks the parameters can be calculated nu- 
merically. In particular, for five-spin blocks we have, in the XX 

limit (for which analytical results are available), 75 = — 
0.577, and in the Heisenberg chain 75 ~ 0.512. 

In blocks with an odd number of spins, the original spins 
Sj (i — l,...,n) are represented by the effective spin So, 
with 'weights' given by the coefficients (? in and cj n , defined 
through the operator identities (valid in H) 



Sf = <„Sq and Sf = c z in S z Q . 



These are useful in the calculation of correlation functions. 
Note that c\ n = c x n n = y„ and c\ n = & n n = 8„y„. For three- 
spin blocks we have 



1 



c 2,3 = 



! + M A °+\/ A o + 



and 



3 = -A A + jAl + 8 )c 



'2,3- 



'»() ' 



Equations JA4> . iA6\ and iA7\ constitute the recursion re- 
lations defining the RG steps in the MDH scheme. 
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Appendix B: SECOND-ORDER CALCULATION OF 
CORRELATION FUNCTIONS 

Let us assume that a two-spin block, as the one shown in 
Fig. |22j is selected for renormalization at some point of the 
RG process. In terms of the states of Si and S2, the eigenstates 
of the block Hamiltonian Ho, with the corresponding energies, 
are 

1*0) = (|U) - lit)) , E Q = - (| + iAo) Jo, 
|^i) = ITT), 1^2) = |U), £i=£ 2 = |Ao7 , 



and 

I^3) = ^(|T4) + I4T)), e 3 = (±-±Ao)Jo. 

The connection between the two-spin block and the rest of the 
chain, through the neighboring spins 5/ and S r , is described 
by the Hamiltonian Hi r in Eq. (lAli . 

Denoting by |A,-) the states of all other spins in the chain, 
and assuming that in the thermodynamic limit there is a 
unique ground state |Ao), the ground state of the whole chain 
can be written, at zeroth-order in perturbation theory, as \go) — 
|Ao, *?())■ Up to second order in J^ r /Jo we obtain a corrected 
state 



\g) 



i MO 
i,j k^O 



(E -E k )(E -E e 



{Aj^k \Hi r \Aj,m ) (Aj,V \Hi r \A ,m ) 
(Eo-E k ) 2 



(Bl) 



A second-order estimate of the expectation value of any oper- 
ator O is simply given by 



_ (g\o\g) 



i\g) 



For the calculation of correlation functions involving spins 
in the block, we write O — O^Oa, where and Oa are 
operators acting on the subspaces defined by the states |*P,-) 
and |A/), respectively. Expanding Eq. iBl\ . we obtain an ex- 
pression for (g \0\g) with terms containing combinations such 
as (Wi \0<ji\ *¥j) and (Aq \S^Oa \ Ao), which is rather cumber- 
some to write here. As examples of the final results obtained 
in the Heisenberg limit, we have 



S n being any spin other than Si, S r , Si and S2. These 
expressions depend explicitly on expectation values like 
(Ao |S/ • S,-| Aq), which on their turn depend on expectation 
values involving spins neighboring the blocks in which Si and 
S r will be decimated. As the RG proceeds, this generates a 
hierarchical structure, which can be solved backwards by as- 
suming that the correlation between the spins in the very last 
block to be renormalized takes it zeroth-order value. It is inter- 
esting to notice that the correlation between two spins which 
are not decimated in the same block is at most of order Ji r jJo- 
Similarly, in the XX limit we have, for instance, 



l\g) 



-1 



= 1 



I 2 
1 J i 



-Jl-4 J -£{A \SfS x r \A ), 



i\8) 



-1 



1 



3 J i 



Jf+J 



16 



Jl 



±^(A |S,-S r |Ao), 



(Si-S 2 ) 



-h 



S/-Si) = 



1 J l 



Jf + Ji 



16 



+ i^(A |SrS r |Ao) 



i + i^)£<*>l**l*>> 

n +4 J ho 



1 /, J rJ r A^l^r 



4^<Ao|S?S;|A )), 



-g (Ao Ao) - ^ (Ao \S&,\ Ao)) 



and 



(S n -Si) g = -\g 



l + ||j|<Ao|S„-S/|Ao) 



J A Jr 



4 Jo) Jo 



(A |S„-S,-|Ao) 



S n S l.2) =8 J Pt( A 



< Jo 



Notice that expressions for the zz correlations may involve 
other expectation values of both xx and zz correlations. 
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